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COMPUTATION OF TURBULENT NEAR WAKE 
FOR ASYMMETRIC AIRFOILS 


George S. Deiwert 
Ames Research Center, NASA 
Moffett Field, California 94035, USA 


A numerical procedure for studying the turbulent near wake of two-dimensional 
airfoil sections is presented. The Reynolds averaged Navier-Stokes equations 
are written for flow about bodies of arbitrary geometry and solved on an 
arbitrary nonuniform curvilinear computational mesh. Eddy viscosity and Rey- 
nolds stress turbulence transport models are considered. Specific examples 
are shown for the RAE 2814 airfoil section by using an algebraic eddy viscos- 
ity model with streamwise relaxation and the interactive Reynolds stress model 
proposed by Bradshaw. 

Symbols 

c chord length 

Cp specific heat at constant pressure 

t unit vector in Cartesian space 

c y specific internal energy 

3 unit outward normal vector 

p pressure 

5 surface area 

T temperature 

t time 

u,v Cartesian velocity components 

vol volume element of fluid 

Y isen tropic exponent 

6 boundary layer thickness 

p coefficient of molecular viscosity 

p mass density 

Subscripts 

BL boundary layer model 

DS evaluated at dividing streamline 

i laminar 
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t* trailing edge 

v wall 

UK wake model 

6 evaluated at boundary layer edge 


Suneracrlpta 

* curvilinear space component 

1 fluctuating component 


Introduction 

The prediction of flow fields over airfoil sections with viscous-inviscid 
interactions has received considerable attention over the past several years. 
One region of particular concern is the trailing edge of the airfoil where 
predictions have been most difficult. In this region, the interaction in the 
near wake of the two turbulent shear layers from the upper and lower surfaces 
with one another is important. Each shear layer experiences a different flow 
history and, in the near-wake region, will sustain normal pressure gradients. 
The turbulent mixing process in the near wake is not yet well understood. 
Simplified numerical analyses initiated just at the trailing edge, not 
accounting for normal pressure gradients, may be Inadequate regardless of the 
degree of sophistication of the turbulence transport model used. 

In this study, the two-dimensional Reynolds averaged compressible Navier- 
Stokes equations are used to describe the entire flow field. In this way, it 
is possible to account for both the upstream history of the shear layers and 
the viscous-inviscid interaction at the trailing edge. Differential turbulent 
transport models can be used to describe the Reynolds stresses in the attached 
shear layers and their carry-over Into the free shear flow of the wake. 


Both eddy viscosity and Reynolds stress turbulence transport models are con- 
sidered. Because the eddy viscosity models relate turbulence transport 
directly to gradients of the mean flow, it is not clear that they can ade- 
quately describe the complex interaction process sustained in the vicinity of 
the near wake. For example, in the near wake, the point of zero velocity 
gradient does not generally coincide with the point of zero shear. Also, it 
is not clear how much fluid dynamic detail is necessary to determine the aero- 
dynamic performance of airfoil sections. Representative models of each type 
are considered in this study: (1) an algebraic eddy viscosity model of the 

type suggested by Smith and Cebeci,* which is modified by a relaxation formula 
in an attempt to account for upstream history effects, and (2) a Reynolds 
stress model of the type suggested by Bradshaw et al. The Reynolds stress 
model is used in conjunction with the interaction hypothesis suggested by 
Bradshaw^ in an attempt to deal directly with the interaction in the near-wake 
region. This model, relatively untested for asymmetric wakes, treats the near 
wake as a complex shear layer composed of two overlapping simple shear layers 
with shear profiles of opposite sign. This model accounts for upstream his- 
tory effects in a natural manner and it is hoped it can support analyses of 
turbulent near wakes. By studying both models and retaining a rigorous set of 
equations and an adequate computational control volume, it may be possible to 
assess the sensitivity of airfoil computations to the details of the near-wake 
interaction and also assess the sensitivity of near-wake computations to the 
sophistication of the turbulence transport model. 


Computations are made with the strong conservative form of the time-dependent 
compressible Navier-Stokes equations by using MacCormack's finite volume mixed 
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explicit/implicit/method-of-characteristics algorithm, 4 The equations are 
written and solved In Cartesian x-y space for computational meshes of arbi- 
trary geometry. An exponentially stretched grid with thickness the order of 
the shear layer thickness is used to resolve the viscous flow regime. Near 
the solid surface* this grid is fine enough to resolve the viscous sublayer 
where velocity varies linearly with distance from the surface. Outside this 
viscous region* a second exponentially stretched grid is used to resolve the 
invlscid field to distances far from the body surface (typically IS chord 
lengths). This kind of grid construction supports all the flow field features 
and also permits computational efficiency. To facilitate adequate resolution 
of the near wake* the computational grid is dynamically remeshed during the 
course of the solution so that the mesh always follows the path of the wake 
(as defined by the locus of minimum velocity) . Solutions are advanced in time 
until a steady state is reached. 

Detailed data is sparse for asymmetric near wakes. In this study, the RAE 
data reported by Cook 5 are compared with our computation. Measurements in- 
clude boundary layer and wake pitot surveys and extracted skin friction and 
integral parameters. Comparisons are made with lift coefficient, pressure 
distribution, and velocity profiles. 

Turbulence Models 

Eddy viscosity models incorporate turbulent transport into the molecular 
transport stress tensor by adding the scalar eddy-viscosity transport coeffi- 
cient e , thereby relating turbulent transport directly to gradients of the 
mean-flow variables. In a Cartesian coordinate system, the two-dimensional 
molecular stress tensor can be written as 

Tjj * (P + 0 x^ e x e x T xy e x e y T yx e y^x (p + °y) e y e y (1) 

where the components are defined by 

a x ■ -2y(3u/3x) - A(3u/3x + 3v/3y) 
o y * -2y(3v/3y) - A(3u/3x + 3v/3y) 

T X y * iy x ■ -y (3u/3y + 3v/3x) 

and 

A - -2y/3 


(5) 

( 6 ) 

(7) 

( 8 ) 


(9) 

Typically, the eddy conductivity coefficient is related to the eddy viscosity 
coefficient via a turbulent Prandtl number Pr t where 

Pr t - p c p e/k t (10) 

The simplest eddy viscosity models are algebraic. In our examples, we used an 
algebraic model of the type suggested by Smith and Cebeci where in the boun- 
dary layer near the solid surface we have 


The total shear (molecular plus turbulent) is written as 

T - + T t - (P + ° x )e x e x + T xy e x e y + T yx e y e * + (p + c y )e y e y 

where 

o x ■ -2(y + pe)(3u/3x) - X(3u/3x + 3v/3y) 
o y « -2(y + pc)(3v/3y) X(3u/3x + 3v/3y) 

t xy - T yx « -(y + pe) (3u/3y + 3v/3x) 

In a similar manner, turbulent heat transport is defined in terms of mean 
energy gradients and an eddy conductivity coefficient k £ so that 

Q t - k t VT 


( 2 ) 

(3) 

(4) 



c •» t 2 l(3u/3y) ? + (lv/8x) 2 J^ 
i » 0.4 nil - #xp(-!n/A)] 

A » 26 )J v //pyT w 


And in th« out«r part of the boundary layer and in the wake we have 


e 


n 


0.0168 u fi 6J/{l + 5.5[(n-nDS)/«l 6 > 

I (1 - u/u fi )dn 

^DS 
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(ID 

( 12 ) 

(13) 

(14) 

(15) 


To approximate the influence of upatream history a simple relaxation procedure 
is used such that 


€ (£, n) ■ ac(5 - 4C. n) + (l - a) e e q(£» n) (16) 

where e eq is defined by equations (11) and (14), a is a relaxation param- 
eter with value between zero and one and is the local streamwise computa- 

tional mesh spacing. For a ■ 0, there is no relaxation and for a ■ 1, the 
eddy viscosity is frozen. A value of 0.3 is used for o in the examples. 

The turbulent Prandtl number is assumed constant at 0.90. 

This model has been used in previous studies 6 * 7 and found to support complex 
interacting flows reasonably well qualitatively. However, it is doubtful that 
a model developed for boundary layers and wakes will describe the complex tur- 
bulent mixing in the near wake of lifting airfoils. 

Reynolds stress models require separate transport equations to describe the 
components of the Reynolds stress tensor. Their dependence on the mean flow 
gradients is indirect as opposed to the direct dependence in eddy coefficient 
models. The simplest Reynolds stress model is the one-equation, thin-shear- 
layer model. In this example, the thin-shear- layer Reynolds stress model pro- 
posed by Bradshaw et al. 2 is considered. The turbulent transport in the near 
wake is described using the interaction hypothesis formulated by Bradshaw et 
al.'* This concept was first applied to symmetric wakes by Morel and Torda 6 
and more recently by Huffman and Ng.° In this approach, the two shear layers 
from the upper and lower surfaces are considered as separate and distinct 
"simple" shear layers. Experiments 10 indicate that, in the near wake, the 
shear layers interact with each other only in a very thin interfacing region. 
For the most part, they continue to behave as a flow in which the shear stress 
does not change sign and which has one set of large eddies only. These two 
simple shear layers then "time share" an ever increasing region as they pro- 
gress downstream from the trailing edge. Figure 1 is a schematic illustration 
of the free interaction concept in the near-wake region. The resultant shear 
stress in the interaction region can be determined by superposing the 
stresses resulting from the two simple shear layers. The attractive feature 
of such a concept is that the complex flow should not be very different from 
those for the corresponding simple shear layers alone. 


The transport equation for Reynolds stress is written in a coordinate system 
corresponding to the thin shear layer as 


3t+ 3ut' i ~ 9vt' t ' 

at as an 


C + T + 


' 3u _ /|t+/p n _ V T f 
.3n ~ L + 


( 17 ) 


where t + - - p u’v’ 

The t erm c*T + 3u/an represents generation (sometimes called production), 

— c + i + /| T + /p f/L + represents destruction (dissipation), and -a-r + V x + /3n is 
transport by bulk diffusion. 
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For certain types of simple shear flows, universal functions have been pro- 
posed for c+, L+ and V T +. Generally, c+ Is assumed constant and equal to 
0.30. Typically, it represents twice the ratio of shear stress t+ to turbu- 
lent kinetic energy. For boundary layer flows Bradshaw et al. 2 suggest the 
functions shown in figure 2 for L* and G where 


V T + - cVItVpWG (18) 

These functions are represented analytically in the present study by the fol- 
lowing relations: 


L+/5 


0.40 n/6 , 

- 0.04 + 0.22(1 -n/6)n/6 , 

816.5 exp(-9.866 n/6) • 


Ga 6 //|x+/p|^ 


(1.117 + 16.94 n/6)n/6 
• 90.87 n/5 - 49.73 
18.69 n/6 + 14.86 


n/6 s 0.1818 
0.1818 < n/6 S 1.10 
n/6 > 1.10 
, n/6 < 0.6286 

, 0.6286 £ n/6 s 0.895 

, n/6 > 0.895 


(19) 


( 20 ) 


In this study, the length scale is modified for small values of n by a van 
Driest 11 damping term so that 


L + /6 - 0.40(n/6) [1 - exp(-n/A)] , n/A < 4.0 


For fully developed wake flows Morel and Torda 8 suggest the functions shown in 
figure 3 for L and G. These functions are represented analytically in this 
study by the following relations 


L + /6 


816.5 exp{-9.866[l-(n/6)] } , 
■0.087 , 

816.5 exp(-9,866 n/6) , 


G 4 W/iTW max 


-20 

34 n/6 - 20 
■ 31 n/6 - 19.3 
20 n/6 - 12 
2 n/6 + 4.4 


n/6 < 0.07289 
0.07289 S n/6 < 0.92711 
n/6 > 0.97211 
n/6 £ 0 

0 < n/6 < 0.2333 
0.2333 S n/6 < 0.6636 
0.6636 < n/6 < 0.9111 
n/6 a 0.9111 


( 21 ) 


( 22 ) 


where A Uj,^ is the maximum velocity difference across the wake. 

In the near wake, the empirical boundary- layer length and velocity functions 
are relaxed to the empirical functions for a fully developed wake using a 
relaxation length that is a multiple of the average of the two boundary layer 
thicknesses at the trailing edge. Hence, we have for the near wake 


L + /6 - #(L + /6) BL + (1 -tfULVd)^ 

(23) 

V T + - «(V t + ) B l + (1 - *)(V t + )wk 

(24) 

9t - exp(- 4 S te /8 6 te ) 

(25) 


where 4 £ te is the distance from the trailing edge, 6 te is the average of 
the upper and lower boundary- layer thicknesses, and 8 is the length scale 
multiplier. Experimental studies on turbulence relaxation length scales sug 
gest a value of 10 for 8, which was used in these studies. 


To describe both shear layers simultaneously, two equations are necessary. 
Letting equation (17) represent the shear layer where the stress is positive 
we can write the stress equation for the shear layer where the stress is 



negative as 


6 


3t" 


9ut' 


3vt‘ 


. ,tE7pT| . HZf 

[3n L" 3n 


(26) 


3t 35 3n 

where c~ • -c+, L”/6"(n/6’) ■ L+/6+(n/6 + ), Vf * c"/|t"/p . 

The Reynolds shear stress is then determined from the sum of the two compon- 
ents so that 

T+ + T- (27) 


A | A | 

- p u'v' 


A I A I 

- p V u 


In regions where there Is no Interaction, e.g., boundary layers and the outer 
region of the wake, one component will be zero. Where there Is an interac- 
tion between the two shear layers, as in the near wake, the superposition 
principle will apply. 

The normal Reynolds stress components are approximated by 

-P 0^ - -5(It+I + !t-|) 

-pO^ - -2(|t+| + |t-|) 


(28) 

(29) 


based on the assumptions that the turbulent kinetic energy is always parti- 
tioned 5:2 and that the ratio of shear to kinetic energy c + /2 is approxi- 
mately 1:7. 

Reynolds analogy is used to relate the turbulent heat flux terms to the Rey- 
nolds stress components so that 

-Cpv’T' - (t+ + x") Y(AejMu) fi (30) 

-CpfiTF - -5 (| t+I + Ix-I) Y(A ei /Au) 6 (31) 

where (Ae*)$ and (Au)$ are the differences between the local values of e^ 
and u ana their values at the outer edge of the shear layer. 

In this initial study, no attempt is made to include extra rates of strain 
resulting from streamline curvature or other causes. These effects can be 
large and may be Important in asymmetric near wakes when there is a strong 
vi8cou8-invi8cid interaction. Still to be studied, as well, is the consider- 
ation of a differential-length scale model. The application of ti.e interac- 
tion hypothesis to near wakes is relatively recent and the inclusion of these 
additional levels of sophistication are perhaps best considered later one at a 
time. 


Solution Procedure 


The equations describing the flow field are written in integral form as 


► 

3t 

vol 


U d vol + J 
S 


H*n ds * ♦vol 


(32) 


where for Cartesian momentum components 


p 


P* 

pu 


puq + T-e x 

pv 


pvq + T*e y 

e 

* H - 

eq + T*q - kVT - Q t 

T + 
t" 


T+q 


U 


(33,34) 
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# » 


0 

0 

0 

0 


c+r+(3u/3n - /|t+/pI/l+) - 3T+v T +/3n 
c"T*(3u/3n - /(t-7p7/L-) - 3r"V T -/3n, 


(35) 


The equations below the dashed line are considered only with the Reynolds 
stress model. They are written in terms of the curvilinear coordinate system 
corresponding to the thin shear layer where g is the 8treamwi.se direction 
and n is the transverse direction. Although the momentum components are 
written in the Cartesian coordinate system, the velocity and heat flux vectors 
are written in the curvilinear mesh coordinate system so that 


q 

VT 


Qt 


ug ? + vg n 

(36) 

(3T/3g)I c + (3T/3n)l n 

(37) 

k t V T , e model 


< 

(3.r, 

(-Cpu’T')g 5 + (“Cp ■O'T’)^ , T t model 



The shear stress is given by equations (5) through (8) for the eddy viscosity 
model. For the Reynolds stress model T ■ t£ - T t where T ^ is defined in 
equations (1) through (4) and 


T t - C-pO^)t 5 i 5 + (-P u'v’)l t i n + (-p« 7 u T )i n J ? + (~pFT)| n J n (391 
and the Reynolds stress components are defined by equations (27) through (29) . 


Equation (32) can be solved for arbitrary geometries with computational meshes 
of arbitrary configuration. For simplicity in treating the Reynolds stress 
equation, the coordinate system should be orthogonal and body oriented in the 
viscous region. This kind of coordinate system also simplifies the character- 
istic operator in the MacCormack** method used to solve the equations. 

Constructing the transformation between Cartesian space and the computational 
space we have 

£ - S(x,y) , n - n(x,y) (40) 

with a Jacobian of the transformation 

J » (3x/3g)(3y/3n) - (3y/3g) <3x/3n) 


The covariant base vectors can be written as 

- (3x/3g)e x + (3y/3g)S y 

► 

ilr, ■ (3x/3n)<£ x + (3y/3n)2y, 
and the contravariant base vectors as 

g£ « (3g/3x)e x + (3g/3y)e y 

► 

gn - On/Sx)^ + On/3y)e y . 

The transformation metrics are 

g ec - " Ox/3g) 2 + Oy/3g) 2 

g nn - g n -g n - (3x/3n) 2 + (3y/3n) 2 
g gn “ g ng “ 0*/H)Ox/3n) + Oy/aO Oy/an), 


(41) 


(42) 


( 43 ) 


and 
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g« . |C.$C . oe/ax) 2 + (35/3y) 2 
g nn . jn.j' a On/3x) 2 + On/3y) 2 
g n£ • g Cn - (3g/3x)On/ax) + (3£/3y) On/3y), 


( 44 ) 


For che computational coordinate system we can write equation (32) in differ- 
ential form as 


3U J_ 3(H*ng A Sg) j_ 8(H»it n A S n ) 

3t l a as Lj an “ 


(45) 


where L* - vol^j/AC , Lj ■ vol^j/An . The flux components can be written as 


p u 

|puu + T *$*•$?> 


pvu + T*iy»a^ 

(H’i^AS^) - JAq eu + - kVT*$£ - Q t *$* 


(46) 


(H-n n 4S n ) 



( 47 ) 


where > 

i^AS^ - (g^/t^gU) (/g nn An) » JAnS^ (48) 

and 

n n AS n - (g n /*^ ?fT 0(* / i^'AC) ■ JAtI n (49) 

For computational efficiency in the viscous region, the term (H*tl T1 AS n ) is 
split into a parabolic part, which is treated implicitly, and a hyperbolic 
part, which is treated explicitly. Thus for the parabolic part we have 


(H*n n A S n ) p - JAC 


(T - p) •2 x *f n 
(T-p)*e *f n 

(T-p)*$*$ n - kVT*$ n - Qt*I n 


( 50 ) 


and the hyperbolic part 


(H*a n A S n ) H - JA£ 


pv 

puv + p?^*^ 
pvv + p£y*$ n 
(e + p)v 


T + V 
T "V 


( 51 ) 


For the Bradshaw Reynolds stress model, there are no gradient diffusion terms 
and the model transport equations are hyperbolic. In practice, the bulk dif- 
fusion velocities V T + and V T - are combined with the convection velocity v. 
The source term ♦ then contains only the generation and destruction terms. 
This source term can lead to numerical instability if sufficient care is not 
exercised in its treatment. Linear stability analysis indicates the equations 


are unconditionally unstable. In the steady state, an exact balance Is main- 
tained between generation, destruction, bulk diffusion, and convection. In 
the transient case, however, there will be regions in the flow where turbu- 
lence Is being generated faster than it is being carried away or dissipated. 
Artificial constraints must be impose i on the magnitude of the generation and 
destruction terms so there is no unstable exponential growth or decay of the 
Reynolds stress. In this study, this is achieved by artlflcally limiting the 
magnitude of generation less destruction such that 

j aa/3n - ✓Ir/p | /Lj < 0.02/At 

where At is the time advanced by the differential operator. Although some 
time accuracy is lost, the growth limitation is not altogether inconsistent 
with the implicit treatment of the other viscous terms. 

In addition to the artificial constraint on the magnitude of <J>, the distribu- 
tion of $ across the shear layer is quadratically smoothed to minimize the 
development of unstable oscillations. To further minimize the chance of in- 
stability, an eddy viscosity transport model is used to obtain starting solu- 
tions, thus avoiding problems associated with a large initial transient 
change . 

Mesh Construction 

A computational mesh is constructed around the airfoil geometry as follows. 
First, an orthogonal body-oriented mesh is constructed around the airfoil 
(with the airfoil surface treated as one mesh line) and extended out a dis- 
tance sufficient to capture the viscous-dominated flow region (typically a few 
percent of the chord length of the airfoil, say 10%). This mesh is extended 
downstream from the trailing edge of the airfoil (typically 15 chord lengths) 
at an angle relative to the mean chord line corresponding to the angle of at- 
tack of the airfoil to the oncoming stream. The mesh spacing normal to the 
airfoil surface is exponentially stretched with the first mesh spacing set 
equal to 1/1 2th of the chord length divided by the chord Reynolds number. 

There are 20-30 mesh points distributed across this viscous layer to the outer 
edge of the zone. Points are generated (or input) along the airfoil surface 
and distributed exponentially downstream in the wake to a maximum total of 100 
"surface" points in the wake and on the airfoil. Second, a larger and coarser 
mesh is constructed extending from the outer edge of the viscous mesh to a 
distance far from the airfoil (typically 15 chord lengths). Again, exponen- 
tial stretching is used away from the surface and some 20 to 30 points are 
distributed across the region. The radial mesh lines are required to be con- 
tinuous with those in the viscous mesh, but the requirement of orthogonality 
is relaxed so that mesh lines neither converge (cross) nor diverge too rapid- 
ly. The mesh constructed for the RAE 2814 airfoil is shown in figure 4. 

As the flow develops, the position of the viscous wake moves with time. To 
assure the capture of the wake in the viscous mesh, the mesh downstream from 
the trailing edge of the airfoil is dynamically adjusted during the course of 
the solution in a direction normal to the mean chord of the airfoil so that 
the slope of the wake mesh matches the flow direction in the core of the wake. 
This dynamic remeshing is initiated sometime after the impulsive start (typi- 
cally after a dimensionless time of 1) to allow time for the starting vortex 
to be shed downstream. Mesh adjustment is realized gradually with the adjust- 
ment never exceeding the smallest transverse mesh spacing in any one time step. 

Results 

Computed results have been obtained for the RAE 2814 airfoil section corres- 
ponding to the test conditions of Cook 5 where the free stream Mach number is 


0.725, the angle of attack is 1.44°, and the chord Reynolds number is 15 * 10* , 
Solutions were generated from an inviscid impulsive start for a dimensionless 
time (U^t/c) of approximately 0.2. The eddy viscosity model was then inclu- 
ded and the no-slip boundary condition imposed, and the solution continued to 
a dimensionless time of 12. Comparisons of computed and experimental wake 
velocity profiles just aft of the trailing edge are shown in figure 5. In the 
first 10% of dimensionless downstream distance, the experimental data show 
fuller profiles from the upper surface than the computations, and less full 
profiles from the lower surface than the computations. Further downstream 
the agreement is quantitatively very good and symmetry is rapidly approached 
in both computation and experiment. 

The computation was extended to a dimensionless time of 13 using the Reynolds 
stress model. Comparisons of computed velocity profiles in the near wake, 
obtained with eddy viscosity and shear stress models, are shown in figure 6. 
Just aft of the trailing edge, the Reynolds stress model solution yields an 
even less full upper layer profile and a fuller lower layer profile than the 
eddy viscosity model solutions. These differences diminish rapidly as the 
flow progresses downstream and, for the upper layer, the trend actually re- 
verses. This suggests too large a Reynolds stress relaxation length in the 
central wake region. A length of 10 6 te was used in the computations. 

Figure 7 shows comparisons of near wake shear profiles for both the eddy vis- 
cosity and Reynolds stress computations. The distinctive difference in the 
two distributions appears in the central core of the wake where the r ddy vis- 
cosity model exhibits large values of sh^ci. stress and a very large gradient 
where the stress changes sign. These large differences in turbulent shear 
distribution and magnitude, however, seem to have little effect on the mean 
velocity profiles, as shown in figure 6. The influence of these differences 
in shear distribution and magnitude is even less on the pressure distribution 
over the airfoil and in the near wake. Figure 8 is a comparison of computed 
and experimental surface-pressure distributions. Except for the vicinity of 
the trailing edge, the agr^ment between computation (with both turbulence 
models) and experiment is very good. There is some streamwise oscillation in 
the computed pressure distribution near the trailing edge and the values are 
typically smaller than the experimental results. The difference in magnitude 
of pressure coefficient is coupled with the velocity distributions at the 
trailing edge where the computed displacement thickness on the upper surface 
is greater than the experimentally determined thickness. Although the near- 
wake pressure prediction is slightly better using the Reynolds stress model, 
there is still room for further refinement before good quantitiatve agreement 
can be reached. 

Concluding Remavxs 

A numeri.’l. scheme has been described to study interactive flows in the vicin- 
ity of two-o.'mensional near wakes. The Reynolds averaged Navier-Ftokes equa- 
tions are used to ensure proper treatment of viscous-invisc id interactions, 
and both eddy /iscosity and Reynolds stress turbulence transport models are 
considered to describe the complex mixing in the viscous-viscous interaction 
region. Altlvugh differences in turbulent shear stress distributions in the 
near wake for the two different models appears to be large, the Influence on 
pressure distribution and resulting lift and drag on the airfoil is small. 

The greatest source of disagreement between computed and experimental pressure 
distribution is related to the large computed displacement thickness just 
ahead of the trailing edge. More detailed turbulence transport models may be 
necessary in the wal 1-shcar-f low region to realize improvement at the trailing 
edge . 
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Fig. 1. Schematic of the over- 
lapping shear layers at asymmetric 
trailing edge. 



Fig. 2. Length and velocity scale 
functions for attached wall layers. 



Fig. 3. Length and velocity scale functions for fully developed wakes. 
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Fig. 4. Computational grid for 
RAE 2814 airfoil section. 
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Fig. 5. Wake velocity profile 
development for RAE Section 2814; 
M = 0.725, a = 1.44, Re = 15 x 10 
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Fig, 7. Computed near-wake turbulent shear profiles; eddy viscosity 
and Reynolds stress models. 
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Fig. 8. Surface and near-wake pressure distributions for RAE Section 
2814; M = 0.725, a = 1.44°, Re = 15 * 10 & . (a) Surface distribution, 
(b) Near-wake distribution. 








